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Abstract 



r-j 

In the flavor basis there are seven cases of two vanishing minors in the neutrino 



mass matrix which can accommodate the present neutrino oscillation data including the 



recent T2K data. It is found that two of these cases, namely B$ and Bq predict near 
maximal atmospheric neutrino mixing in the limit of large effective neutrino mass. This 
feature remains irrespective of the values of solar and reactor mixing angles. A non-zero 
■ reactor mixing angle is naturally accommodated in these textures. 

o 
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1 Introduction 

During the past decade there has been considerable experimental development in the de- 
termination of neutrino masses and mixings pp. Recently, T2K experiment [2] has given 
unambiguous hints of a relatively large 1-3 mixing angle. In this light, it is natural to look 
for models which, naturally, accommodate a non-zero value of reactor mixing angle while the 
atmospheric mixing angle remains near its maximal value. Recently many papers have ap- 
peared which reproduce the relatively large value of the reactor mixing angle [3]. 
There are mainly two approaches to explain neutrino mixings: 1) Mass independent textures 
[I] which lead to mixing matrices independent of the eigenvalues. The most celebrated ex- 
ample of this category is the tribimaximal (TBM) [5] scenario which has been derived from 
family symmetries and predicts a vanishing 1-3 mixing angle #13 = 0, maximal 2-3 mixing 
angle #23 = vr/4 and 1-2 mixing angle 612 = sin _1 (l/\/3). Non-zero #13 can be accommodated 
in TBM and other similar models by considering deviations from symmetry. 
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Class 


Zero Minors 


Ai 


C3,3> Cs,2 


A 2 


C 2 ,2, Cs : 2 


B 3 


C3,3) Cs,l 


B 4 


C2,2, C2,l 


B 5 


C3,3) Cl,2 


B 6 


£-2,2) Cl,3 


D 


C3,3) C*2,2 



Table 1: Experimentally allowed classes of two zero minors, here C# denotes the zero minor 
corresponding to the element of M v . 

2) Mass dependent textures which induce relations between mixing matrix elements and mass 
eigenvalues. Such textures naturally accommodate a non-zero #13. Some examples of these 
are zero textures [5], vanishing minors [7J [5], hybrid textures [9]. Zero textures have been 
particularly successful in explaining both the quark and the lepton masses and mixings. 
In this work we identify a class of mass dependent textures which supplemented with the 
assumption of a large value of effective neutrino mass M ee naturally predict near maximal #23 
and non-zero #13. Recently, it was shown by Grimus et al [10J that near maximal atmospheric 
mixing is predicted for class B3 and B4 of two zero textures supplemented with the assump- 
tion of quasi degeneracy. We consider two vanishing minors of the neutrino mass matrix in 
the flavor basis together with the assumption of large M ee . This assumption is well motivated 
by the extensive search for this parameter in the ongoing experiments. We found that the 
two cases of two vanishing minors viz. B§ and Bq [Table 1] in the limit of large M ee predict 
near maximal atmospheric mixing and this property holds irrespective of the values of solar 
and reactor mixing angles. The seesaw mechanism [TTJ is regarded as the prime candidate for 
understanding the scale of neutrino masses. In the framework of type-I seesaw mechanism, 
the effective Majorana mass matrix M v is given by 

M v = -MdM^MI (1) 

where Mo is the Dirac neutrino mass matrix and Mr is the right-handed Majorana mass 
matrix. In the framework of type-I seesaw mechanism M u is a quantity derived from Md and 
Mr. Therefore, zeros of Mr, and Mr have a deeper theoretical meaning. In a basis where Mb 
is diagonal, the zeros of Mr propagate as zero minors in M v . 

2 Symmetry realization 

In the basis where the charged lepton mass matrix is diagonal, there are fifteen possible two 
vanishing minors in M v . Out of these fifteen only seven patterns [Table-I] viz. Ai, A 2 , B 3 , 



2 



£?4, B 5 , Bq and D can accommodate the neutrino oscillation data. Of all the allowed two 
zero minors in the neutrino mass matrix only three cases B 5 , B e and D provide non-trivial 
zero minors, all other cases reduce to two zero textures when confronted with the neutrino 
oscillation data. We work in a basis where Mp is diagonal (Md = diag(x,y,z)), and the 
neutrino mixing arises solely from Mr. In this basis, a zero entry in Mr propagates as a 
vanishing minor in the effective neutrino mass matrix M v . Here, we focus on B 5 and B§ class 
of vanishing minors. To obtain the classes -B5 and Bq of neutrino mass matrices, we extend 
the Standard Model (SM) by adding three right-handed neutrino singlets v R . and one scalar 
singlet x- m order to enable the seesaw mechanism for suppressing the neutrino masses Mr 
must have the following structures for B 5 and B e : 



M R (B 5 



( a c \ 
Ode 
c e 



M R (B 6 



(a 

b 



f J 



(2) 



leading to the following effective neutrino mass matrices through the seesaw mechanism 

c( .ri/ cclxz \ ^ / < -.1- /»/.!•// —bexz 

aeyz , M V {B & ) = -~ r — ^ 

2 Vf + ae 2 

—adz ) \ 



M V (B 5 ) = 



1 



c 2 d- 



—cexy 
cdxz 



1 1 

c y 



aeyz 



bfxy 
—bexz 



bfxy 
~afy 2 
aeyz 



aeyz 
b 2 z 2 



(3) 

A general procedure for enforcing zero textures in arbitrary entries of the fermion mass 
matrices using abelian family symmetries has been outlined in [T2j • The symmetry realization 
of all the allowed one zero and two zero textures was recently presented in [13]. For the 
symmetry realization of B 5 and B e textures of two zero minors we consider a small cyclic 
group Z3 which corresponds to the minimal group since Z2 leads to a non-diagonal charged 
lepton and Dirac neutrino mass matrix. Under Z 3 the SM Higgs doublet remains invariant 
and the leptonic fields are assumed to transform as: 



Dli ~ > Dl 1: Iri — > Iri, vr x — > vr 1 , 

D L2 -)■ UJD L2 , Ir 2 -> ul R2 , Vr 2 -)■ U)Vr 2 , 

D Ls uj 2 D Lz , l R3 u 2 l R3 , vr 3 -> u 2 v R3 , 



(4) 



where to = e l27r ^ 3 . Hence the bilinears D^.l^ and D^.v^, relevant for Mi and Mp transform 

as 



DidRk ~ D Li u Rk ~ 



/ 


1 




LO 2 


\ 






1 


CO 




V 




U 2 


1 


J 



(5) 



The SM Higgs doublet remains invariant under Z 3 leading to diagonal Mi and M^. The 
bilinear VR.VR k relevant for Mr transforms as 



VR i VR k ~ 



/ 


1 


UJ 


u 2 


\ 




LO 


UJ 2 


1 




V 


LO 2 


1 


LO 


J 



(6) 
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We assume a scalar singlet x transforming as \ — > uj\ for class B 5 and x ~^ ^ 2 X fo r class B e 
which leads to the following Z 3 invariant Yukawa Lagrangians for classes B 5 and B 6 : 

- C {Bb) = Y^Dl^Ir, + Y} 2 D L2 <pi R2 + Yi 3 D L3 ct>l R3 + Y°D Ll 4>v Rl + Y£D L2 4>vr 2 + 

yM y m M M M M 

Y£D L Jv R3 + -^-u T Rl C- x v R3 X + -f^C" W + -^Vr.C- 1 ^ + -^vl 2 C- x v Rz + H.c. (7) 

- £ (Be) = Y^Dl^Ir, + Y^Dl^Ir, + Yi 3 D L3 ct>lR 3 + Y^D L ^v Rx + YgD L ^v R2 + 

yM yM tutM j^M 

Y 33 D L 3 WR 3 + — ^i?^ ^R 2 X + — Vr^ ^flsX + -g— ^Jij + —^—^2^ ^ 3 + H.C. (8) 

where = ir 2 0*. Next, we show how a large effective neutrino mass can arise in such a model. 
We note that M R contains two types of mass terms viz. 1) Bare mass term which does not 
need a scalar singlet and is invariant by itself. 2) Terms arising from Yukawa couplings to x- 
The scale of latter is restricted by the scale of Z 3 breaking while there is no such restriction on 
the bare mass term which can have a higher mass scale. It can be seen from eqn.(3) that the 
ee and fir entries of M v have contributions to their numerators from ee and \xr entries of M R 
which arise from the bare mass term. We assume the mass eigenvalues of Mp to have same 
order of magnitude which leads to a large value of ee and fir entries of M v while the other 
elements of M v are suppressed, thus, leading to a large value of M ee . Since these textures are 
realized at the seesaw scale, the Renormalization Group (RG) evolution of the parameters of 
M v from the seesaw scale to the electroweak scale needs to be taken into account. It is well 
known that the RG effects are most prominent for the quasidegenerate mass spectrum which 
is precisely the case here due to the assumption of large M ee . However, it is also known that 
zero minors in M u , at a given energy scale, remain zero at any other energy scale at the one 
loop level [TJ. This is because the matrices at any two energy scales \i\ and //2 are related 
by M u (fii) = IM v (pL2)I, where / is diagonal, positive and non singular. The operation of 
diagonal matrices from left and right on M v does not alter the zero minors of M v leading to 
zero minors in M v at any other scale. 



3 Formalism 

We reconstruct the neutrino mass matrix in the flavor basis assuming neutrinos to be Majorana 
particles. In this basis, a complex symmetric neutrino mass matrix can be diagonalized by a 
unitary matrix V as 

M u = VMf ag V T (9) 

where M^ %ag = diag(mi, m^, tus). The matrix M v can be parameterized in terms of three 
neutrino masses (mi,m 2 ,m 3 ), three neutrino mixing angles (#12, #23, #13) and the Dirac-type 
CP- violating phase 5. The two additional phases a and j3 appear if neutrinos are Majorana 
particles. The matrix 

V = UP (10) 
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where 



U 



C12C13 



S12C13 



s 13 e- iS \ 



-S12C23 - Ci2S2 3 Sl 3 e J<5 C12C23 - S 12 S 23 S 13 e lS S23C13 



\ S12S23 - Ci 2 C 2 3Si 3 e 



IS 



(11) 



-C12S23 - s 12 c 23 s 13 e ld C23C13 / 



with Sij = sin 8ij and Cjj = cos 6^- and 



P 



I 1 \ 

e ia 
v j 

is the diagonal phase matrix with the two Majorana-type CP- violating phases a, (3 and 
Dirac-type CP-violating phase 5. The matrix V is called the neutrino mixing matrix or the 
Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [15j . Using Eq. (9) and Eq. (10), the 
neutrino mass matrix can be written as 



M, = UPM d ™ 9 P T U T . 



(12) 



The CP violation in neutrino oscillation experiments can be described through a rephasing 
invariant quantity Jcp [16J with Jcp = Im{U e iU^ 2 Ul 2 U^i). In the above parametrization, 
J cp is given by 

Jcp = s 12 s 23 s 13 c 12 c 23 c 2 13 sin 5 . (13) 
The simultaneous existence of two vanishing minors in the neutrino mass matrix implies 



M v (pq)M v (rs) — M v (t u )M u (y W ) = , 

These two conditions yield two complex equations viz. 

3 

(VpiVqiVrkVgk - V t iV u iV vk V wk )mim k = , 

l,k=l 

3 

(Vp'iVq'iVr'kVs'k - V t qV u qV v , k V w > k )mim k = . 



(14) 
(15) 



(16) 
(17) 



Lk=l 



The above equations can be rewritten as 



mi m 2 A 3 e 2ia + m 2 m 3 A x e 2l ^ +8) + m 3 m x A 2 e 2l ^ &) = . 

Ha+W) + m 3mi B 2 e 2i ^ = 



where 



m\m 2 B 3 e ia + m 2 m 3 B\e 



Ah — (UpiUgiUrkUsk — UtiU u iU V kU W k) + (I ■<->■ k) , 
Bh = (Up'iU q 'iU r 'kU s 'k — UfiU u >iU v ' k U w >k) + (I ■<->■ k) , 



(18) 
(19) 

(20) 
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with (h,l,k) as the cyclic permutation of (1,2,3). These two complex eqns.(18) and (19) 
involve nine physical parameters m 1; m 2 , m 3 , 12 , 9 23 , 9%3 and three CP- violating phases a, 
(3 and 5. The masses m 2 and m 3 can be calculated from the mass-squared differences Am 2 2 
and |Am 23 | using the relations 



m 2 = \f mj + Am 2 12 , m 3 = yjm% + |Am| 3 | (21) 

where m 2 > m 3 for Inverted Spectrum (IS) and m 2 < m 3 for Normal Spectrum (NS). Using the 
experimental inputs of the two mass-squared differences and the three mixing angles we can 
constrain the other parameters. Thus, in the two complex eqns.(18) and (19) we are left with 
four unknown parameters mi, a, (3 and S which are, obviously, correlated. Simultaneously 
solving eqns.(18) and (19) for the two mass ratios, we obtain 

mi e - 2ia = A3B1 - A1B3 
m 2 A 2 B 3 - A 3 B 2 



and 



m 3 A 3 B 2 -A 2 B 3 



The mass ratios for class B§ to first order in si 3 are given by 



m 2 ci 2 c^ 3 s 12 

and 



e W ~ 1 _ (25) 



m 3 c 23 4 3 s 12 



The mass ratios for class Bq to first order in si 3 are 



mi e _ 2ia w l _ suc 23 [c 23 e^+s 23 e- , ^ 
m 2 Ci 2 s 23 Si 2 

and 

mi e _ 2i , ^ _c|e^ + c 12 s 13 4 3 (cj^ + e 2 ^ ^ ^ 

m 3 s 23 S 23 Si 2 

In the case of zero textures there exists a permutation symmetry between different patterns 
[T7] . Similarly, there exists a permutation symmetry between patterns B§ and Bq of two zero 
minors which corresponds to the permutation in the 2-3 rows and 2-3 columns of M v . The 
corresponding permutation matrix is given by 



P: 



23 



(28) 



/ 1 \ 

001 . 
V° 1 0/ 

The right-handed Majorana mass matrix Mr for class Bq can be obtained from Mr for class 
B 5 by the transformation 

LR -P 23 M| 5 P 2 T 3 (29) 
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which after the seesaw gives 

= P 23 M*P£ . 

This leads to the following relations between the parameters: 

71 



nB 6 _ nB 5 
u 12 — u 12 i u 13 



13 ' u 23 



--eg*, 5 b ° = 8 b s-tt. 

2 23 



The magnitude of the two mass ratios in eqns.(22, 23), is given by 

P 



m 3 



(30) 
(31) 

(32) 
(33) 

(34) 

2" J \A 3 B 2 -A 2 B 3 ~ J' (35) 
Since, Am^ and |Am| 3 | are known experimentally, the values of mass ratios (p, a) from 
eqns.(32) and (33) can be used to calculate m\. This can be done by inverting eqns.(21) to 
obtain the two values of m\ viz. 



a 



^l e -2ia 

m 2 



while the CP- violating Majorana phases a and /3 are given by 



a 



_\ ar ( A 3 B X -A X B 3 
2 ar9 \A 2 B 3 -A 3 B 2 

1 (A 2 B 1 -A 1 B 2 2iS 
--arg ( -^-^ ^^e zu 



mi = a\ 



1-a 2 



mi = p\ 



Am?, + I Am, 



12 

1 



231 



(36) 



4 Numerical Analysis 

The experimental constraints on neutrino parameters at 1, 2 and 3a [T8] are given below: 

A m 2 _ ? , 8 (+0.22,+0.41,+0.60) ^-5^2 
affl 12 — '• oo (-0. 26, -0.42, -0.59) X iU e V ' 

\l±m 23 \ — ^-<JO(_o.09,-0.18,-0.29) X iU eV ) 

• 2 n _ n o 19 (+0.017,+0.035,+0.052) 
bill t7i2 — U.OlZ(_ 0016 _ 0032) _ 047 ), 

c in 2 /9 - fl /l9(+0.08,+0.18,+0.22) 
Sin tf 23 - U - 4Z (-0.03,-0.06,-0.08)' 

•2/1 _ n not -(+0.007,+0.016,+0.025) / QV \ 
Sin 6»i 3 - U.U^O ( _ 007i _ 0013 _ 020) . \6() 

The observation of neutrinoless double beta (NDB) decay would signal lepton number violation 
and imply Majorana nature of neutrinos, for recent reviews see [19], 120] . The effective Majorana 
mass of the electron neutrino M ee which determines the rate of NDB decay is given by 

M ee = |m lC 2 2 c 2 3 + m 2 s 2 12 c 2 13 e 2ia + m 3 s 2 13 e 2ip \ . (38) 
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Figure 1: Correlation plots for class B 5 , plots (a), (6) correspond to Normal Spectrum (NS) 
and plots (c), (d) correspond to Inverted Spectrum (IS). 

Part of the Heidelberg-Moscow collaboration claimed a signal in NDB decay corresponding 
to M ee = (0.11 — 0.56)e\^ at 95% C. L. [21]. This claim was subsequently criticized in [22] . 
The results reported in [2T] will be checked in the currently running and forthcoming NDB 
experiments. There are large number of projects such as CUORICINO[23], CUORE |24j . 
GERDA [25], MAJORANA [26], SuperNEMO [27], EXO [28], GENIUS [29] which aim to 
achieve a sensitivity upto O.OleV for M ee . In the present work, we take the upper limit of 
M ee to be 0.5 eV [20]. We vary the oscillation parameters within their known experimental 
ranges. However, the Dirac-type CP-violating phase 5 is varied within its full range. The two 
values of mi obtained from the mass ratios p and a, respectively must be equal to within the 
errors of the oscillation parameters for the simultaneous existence of two vanishing vanishing 
minors in M v . The first step in the numerical analysis uses the information of the two known 
mass squared differences along with the constraint of two zero minors and large M ee to get 
predictions for the mixing angles. It is found that both the classes B 5 and Bq predict a near 
maximal atmospheric mixing angle while the other two mixing angles remain unconstrained. 
The atmospheric mixing angle #23 moves towards 7r/4 with increasing M ee as seen in Fig.(l) 
for class B§ and Fig. (2) for class Bq. Thus classes S5 and Bq of two vanishing minors in M v 
naturally predict a near maximal atmospheric mixing angle. 

The second step takes into account the experimental input on the three mixing angles including 
the recent T2K results on the reactor mixing angle. The results for #23 are plotted in Fig. (3) 
for class B 5 and Fig. (4) for class B e . Due to the relatively large value of 9\ 3 , the Dirac-type CP 
violating phase is almost fixed near 7r/2 or 3ir/2 predicting almost maximal CP violation for 
these textures. Fig. (5) and Fig. (6) show the correlation between the two Majorana-type phases 
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Figure 2: Correlation plots for class B e , plots (a), (6) correspond to NS and plots (c), (d) 
correspond IS. 
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Figure 3: Correlation plots for class £> 5 , plots (a), (b) correspond to NS and plots (c), (d) 
correspond IS. 
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Figure 4: Correlation plots for class B 6 , plots (a), (6) correspond to NS and Plots (c), (d) 
correspond IS. 

and M ee : the phases a and (3 approach zero with increasing M ee . As an example, we write 
the numerically estimated mass matrices for the pattern B^, the matrices are obtained for the 
best fit values of Am^ 2 , |Am| 3 |, 6> 12 , #13 given in eqn.(37). For NS we have M ee = 0.11296^, 
9 23 = 43.902°, 5 = 269.463° and for IS we have M ee = 0.1272 e y, 9 23 = 46.058°, S = 89.853° 

0.112956 - 0.000905i 0.000205 - 0.002058i 0.000018 - 0.000166i 
M* r >{NS)=\ 0.000205 - 0.002058i -0.000037- 0.000008i -0.117485 + 0.002109i |, (30.) 
0.000018 -0.000166i -0.117485 + 0.002109i -0.009482 + 0.000114i 

/ 0.127210- 0.001176i 0.000147 - 0.002268i -0.000012 + 0.000177i 
M^(IS)=\ 0.000147- 0.002268i -0.000040 - 0.000006i -0.122682 + 0.002478i |. (4(1) 
\ -0.000012 + 0.000177i -0.122682 + 0.002478i 0.009575 - 0.000154i 

The numerical matrices for pattern B 6 can be obtained from above matrices with the operation 

of 2-3 permutation symmetry. The assumption of large M ee is testable in the ongoing and 

forthcoming experiments [231 EH GHJ EHJ E7J |28j |29] for NDB decay which will either confirm 

or rule out large M ee in the next few years. 



5 Summary 

The recent results of the T2K experiment suggest a relatively large reactor mixing angle. 
Therefore, it is important to look for models naturally accommodating a non-zero value of 
reactor mixing angle while keeping the atmospheric mixing angle near maximal. In the present 
work, we studied the implications of class B§ and Bq of two zero minors in M v for large ef- 
fective neutrino mass. In the context of type-I seesaw mechanism, taking Mi and Mjj to be 
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Figure 5: Correlation plots of Majorana Phases with M ee for class B5 (IS). 
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Figure 6: Correlation plots of Majorana Phases with M ee for class B 6 (IS). 



diagonal, the zeros of Mr propagate as zero minors of M v and the origin of neutrino mixing 
is solely from Mr. We presented the symmetry realization of these patterns using a cyclic 
group Z%. It was found that class B 5 and B 6 predict a near maximal atmospheric mixing 
angle in the limit of large M ee . Furthermore, this prediction is independent of the values of 
the reactor and the solar mixing angles. The assumption of large M ee is testable in the on- 
going experiments for NDB decay since the rate of NDB decay is proportional to M ee . These 
experiments will either confirm or rule out a large value of M ee in the next few years. The 
atmospheric mixing angle approaches 7r/4 with the increasing value of M ee . A reactor mixing 
angle equal to zero is not allowed in these textures, thus, naturally accommodating a non-zero 
#13 as suggested by the recent results of the T2K experiment. Due to the relatively large value 
of #13 the Dirac-type CP violating phase is fixed near 7r/2 or 3tt/2 predicting almost maximal 
CP violation for these textures. 
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